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The cytokine tumor necrosis factor-alpha (TNFa) is one of the key factors during the priming 
phase of liver regeneration as well as in hepatocarcinogenesis.TNFa activates the nuclear 
factor K-light-chain-enhancer of activated B cells (NFkB) signaling pathway and contributes 
to the conversion of quiescent hepatocytes to activated hepatocytes that are able to prolif- 
erate in response to growth factor stimulation. Different mathematical models have been 
previously established for TNFa/NFKB signaling in the context of tumor cells. Combining 
these mathematical models with time-resolved measurements of expression and phos- 
phorylation of TNFa/NFKB pathway constituents in primary mouse hepatocytes revealed 
that an additional phosphorylation step of the NFkB isoform p65 has to be considered in the 
mathematical model in order to sufficiently describe the dynamics of pathway activation in 
the primary cells. Also, we addressed the role of basal protein turnover by experimentally 
measuring the degradation rate of pivotal players in the absence of TNFa and including this 
information in the model. To elucidate the impact of variations in the protein degradation 
rates on TNFa/NFKB signaling on the overall dynamic behavior we used global sensitiv- 
ity analysis that accounts for parameter uncertainties and showed that degradation and 
translation of p65 had a major impact on the amplitude and the integral of p65 phospho- 
rylation. Finally, our mathematical model of TNFa/NFKB signaling was able to predict the 
time-course of the complex formation of p65 and of the inhibitor of NFkB (IkB) in primary 
mouse hepatocytes, which was experimentally verified. Hence, we here present a math- 
ematical model for TNFa/NFKB signaling in primary mouse hepatocytes that provides an 
important basis to quantitatively disentangle the complex interplay of multiple factors in 
liver regeneration and tumorigenesis. 
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1. INTRODUCTION 

Liver regeneration after tissue damage is a tightly regulated spatio- 
temporal process, which is primarily controlled by specific growth 
factors and cytokines (Michalopoulos, 2010). It predominantly 
relies on the fast initiation of hepatocyte proliferation after injury 
and only few proliferative cycles are necessary to restore the liver 
mass and function (Papp et al., 2009). 

Tumor necrosis factor-alpha (TNFa), which is predominantly 
secreted from liver- resident macrophages (Kupffer cells), repre- 
sents one of the earliest stimuli for proper initiation of hepatocytic 
proliferation (Yamada et al., 1997). The most extensively charac- 
terized response activated by TNFa is the highly conserved nuclear 
factor of K-light-chain enhancer of activated B cells (NFkB) path- 
way, which can initiate hepatocytic mitotic cycles after entering 
the nucleus. After partial hepatectomy (PHx), which is the best 
experimental model for studying liver regeneration, the level of 
TNFa increases dramatically within 1 h, leading to fast activation 



of NFkB signaling (Diehl and Rai, 1996). The pivotal relevance 
of TNFa-induced signaling is supported by the fact that in type I 
TNFa receptor- deficient mice, hepatocyte proliferation is strongly 
reduced and liver regeneration is impaired after PHx (Yamada 
et al, 1997). 

Activation of the NFkB axis involves the formation of homo- or 
hetero-dimeric transcription factors composed of different sub- 
units: p50 (NFkBI), p52 (NFkB2), p65 (RelA), RelB, and c-Rel 
(Hayden and Ghosh, 2012). In absence of pathway activation, 
these dimers are disabled by a family of inhibitory proteins includ- 
ing iKBa, IkB|3, IkBs, and iKBy. For example, binding of iKBa to 
hetero-dimeric p65:p50 results in steady- state cytoplasmic com- 
plex retention and therefore prevents dimers from binding their 
respective DNA recognition sites. Upon TNFa stimulation, dif- 
ferent regulatory proteins are recruited in close proximity to 
the TNFa receptor and form the receptor- associated signalo- 
some. Subsequently, IKBa phosphorylation is mediated by a kinase 
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complex composed of IkB kinase-a (IKKa), IkB kinase-^ (IKKp), 
and IkB kinase-y (IKKy or NEMO), followed by proteasomal 
degradation of inhibitory IkBoi as well as release and nuclear 
translocation of p65:p50 complexes. These nuclear heterodimers 
induce transcription of genes involved in the initiation of prolif- 
eration as well as negative regulation of apoptosis (Hayden and 
Ghosh, 2012). 

The modulation of NFkB signaling results in different expres- 
sion signatures in vivo and in vitro (Ashall et al., 2009; Li et al., 
2009; Sung et al, 2009). This could be achieved by posttransla- 
tional modifications such as alternative phosphorylation of the 
involved signaling components. Multiple phosphorylation sites 
have been described for the NFkB subunit p65 (Viatour et al, 
2005). Most common is the phosphorylation of p65 at serine 536 
that is mediated by IKK-|3, Akt, and RSK1 in response to TNFa 
(Sakurai et al, 1999), and ILip (Madrid et al, 2001) stimulation 
and DNA damage (Bohuslav et al, 2004). Other phosphoryla- 
tion sites and kinases are gaining importance for regulating p65 
transcriptional activity. 

The protein kinase C, zeta (PKG;) is capable to phosphorylate 
free p65 at serine 311 when translocated into the nucleus follow- 
ing TNFa stimulation (Duran et al., 2003). Furthermore another 
specific p65 phospho-site, serine 276, has been shown to be a tar- 
get for protein kinase A, c (PKAc), and mitogen -activated protein 
kinase 14 (MAPK14/MSK1; Vermeulen et al, 2003; Jamaluddin 
et al., 2007). These kinases are of particular interest since they 
are localized in the nucleus and provide the possibility to modify 
the transcriptional activity of p65 that translocates to the nucleus 
subsequent to phosphorylation at serine 536 by IKK-|3. However, 
currently there is not much known regarding their direct involve- 
ment in the regulation of canonical NFkB signaling (Joo and 
Jetten, 2008) and the efficiency of antibodies used for detection 
are controversially discussed (Spooren et al., 2010). 

Another important feature of NFkB signaling is ligand- 
independent, basal, and ligand- dependent turnover of the signal- 
ing components. Basal degradation of IkBcx has been reported 
(Krappmann et al., 1996; Pando and Verma, 2000) suggesting 
that two different IkBcx pools, p65-bound and free unbound, 
exist. In both cases, IKK- 13 -mediated IkBcx phosphorylation was 
not necessary for basal protein degradation. However, the pool 
of free IkBoi represents only 15% of the total protein amount 
(Rice, 1993) and thus specific mechanisms are required to explain 
basal degradation of free and p65-bound IkBoi (Mathes et al., 
2008). Combining genetic tools and mathematical models Mathes 
et al showed that free IkBcx degradation is independent of its 
phosphorylation or ubiquitination but mainly mediated by a spe- 
cific C-terminal sequence (PEST) able to fine-tune turnover of 
free IkBoi protein. On the other hand, different from the rather 
instable free iKBa, p65 is relatively stable and upon complex 
formation with IkBcx protects IKBa from fast degradation. IKK- 
P-mediated phosphorylation of p65-bound IkBcx in response to 
TNFa stimulation leads to the release of the complex unmask- 
ing the PEST sequence and results in an efficient recognition 
and fast degradation by the proteasome. Thus, basal and ligand- 
dependent turnover operate on different time scales. Their specific 
contribution to signal amplitude and duration remains to be 
elucidated. 



In the past, the dynamic behavior of NFkB has been intensively 
analyzed in different cancer cells and immortalized fibroblasts 
(Cheong et al., 2006; Ashall et al, 2009). Particularly, time-resolved 
phosphorylation analysis of pathway constituents (including p65 
and IKBa) followed by the activation of target genes involved in 
feedback regulation (e.g., IKBa, A20) leads to the identification of 
sustained oscillatory behavior in NFkB signaling (Nelson et al., 
2004). 

In this regard, NFkB signaling has been extensively studied by 
computational modeling. Carlotti et al. (1999, 2000) presented 
a first computational study focused on the investigation of the 
NFkB shuttling, whereas Hoffmann et al. (2002) studied the whole 
signaling module. The model by Nelson et al. (2004) pointed to 
the importance of the IKBa transcription rate for the frequency 
of oscillations and further analysis suggested IKBa and A20 have 
a major impact on the pathway dynamics (Werner et al., 2008). 
Another study reported a key role for IKBa and IKKs in oscilla- 
tions and protein turnover (Ihekwaba et al., 2004; Cheong et al., 
2006). Recently, the influence of model parameters on the path- 
way dynamics has been studied identifying NFkB concentration 
to be the most decisive parameter for oscillations (Wang et al., 
2011). 

Only few studies included protein synthesis and degradation 
rates as relevant parameters for NFkB signaling dynamics (O'Dea 
et al, 2007). In some cases, IkB turnover is part of the model, but 
NFkB, e.g., p65 turnover has not been included so far, which is 
probably due to the relatively stable levels of NFkB. However, it 
is problematic to predict the influence of stable proteins on sys- 
tems behavior. O'Dea et al. (2007) measured and modeled the 
influence of NFkB binding on the degradation rate of IkB species, 
which was important for their computational predictions. Thus, 
NFkB seems to protect IkB species from unstimulated degrada- 
tion. None of the described studies analyzed and modeled NFkB 
signaling in primary and non-malignantly transformed hepato- 
cytes representing the major cell population of the liver involved 
in regeneration and tumorigenesis. 

Here we present a hepatocyte-specific model for TNFa-induced 
NFkB signaling that considers basal protein turnover and facili- 
tates the prediction of TNFa-induced complex formation of p65 
and IKBa. Since not all parameters of the model are fully identifi- 
able we study an ensemble of 30 different parameter sets that each 
fit our experimental data. Model validation and conclusions are 
based on common properties of the ensemble of models. 

2. MATERIALS AND METHODS 

2.1. EXPERIMENTAL: ANIMAL WORK, ISOLATION OF PRIMARY 
MURINE HEPATOCYTES, AND CULTURE CONDITIONS 

C57B1/6 mice were obtained from Charles River (Wilmington, MA, 
USA) and housed under standard conditions in the DKFZ animal 
facility. All experiments were performed with the approval of the 
German Regional Council of Baden- Wuerttemberg (Karlsruhe, 
Germany) and in accordance with the institutional regulations. 
Isolation of hepatocytes from mice has previously been described 
(Castoldi et al., 20 1 1 ). In brief, six to twelve weeks old animals were 
anesthetized using 10% ketamine hydrochloride (5/100 mg body 
weight) and 2% xylazine hydrochloride (1/100 mg body weight), 
and perfused with HANKS solution supplemented with 0.3 mg/ml 
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collagenase CLSII and 5mM CaCl2. After liver removal, hepa- 
tocytes were plated on collagen-coated culture dishes (1 million 
cells/6 cm BD; Biocoat; Horsham, PA, USA) and cultured with 
adhesion medium (Williams' medium E (Biochrom, Berlin, Ger- 
many), 10% FCS, 0.1% insulin, 100 nM dexamethasone, 2mM 
L-glutamine, and 1% penicillin/streptomycin) at 37°C for 4h. 
After attachment, cells were washed with PBS (GIBCO Life Tech- 
nologies, Darmstadt, Germany) and incubated over-night with 
pre -starvation medium (Williams' medium E, 100 nM dexametha- 
sone, 2 mM L-glutamine, and 1% penicillin/streptomycin). Before 
TNFa stimulation, cells were incubated in starvation medium 
(Williams medium E and 1% penicillin/streptomycin) for 5h. 
Dynamic stimulation of NFkB signaling was achieved by admin- 
istration of recombinant murine TNFa (lOng/ml; R&D Systems, 
Minneapolis, MN, USA) in starvation medium. At different time- 
points after TNFa treatment (5, 10, 20, 40, 60, 120, 180, 240 min), 
medium has been removed and cells were washed with PBS 
immediately before mRNA and protein isolation. 

2. 1. 1. mRNA preparation and real-time PCR 

For isolation of total mRNA from hepatocytes the NucleoSpin 
RNA II kit has been used (Macherey-Nagel Duren, Germany). 
cDNA synthesis (RevertAid H minus, Fermentas, St. Leon-Rot, 
Germany) and semiquantitative real-time PCR (Abgene/Thermo 
Fisher, Epsom, UK) have been performed using SYBR Green ROX 
Mix (Thermo Scientific, Ulm, Germany) according to the manu- 
facturers' protocol. The following primers were used in this study: 
IkBa forward: 5 f -CCT GGC CAT CGT GGA GCA CT-3', IkBa 
reverse: 5'-AGT AGC CTT GGT AGG TTA CC-3'; tubulin for- 
ward: 5'-TCA CTG TGC CTG AAC TTA CC-3'; tubulin reverse: 
5'-GGA ACA TAG CCG TAA ACT GC-3' (Thermo Scientific, Ulm, 
Germany). Cycling program: 95°C/15min followed by 40 cycles 
95°C/15s and 60°C/lmin. Melting curve: 95°C/15s, 60°C/30s, 
95°C/15s. 

2.1.2. Protein preparation, immunoprecipitation f and 
immunoblotting 

Total protein fractions were collected at the indicated time-points 
using the Cell Lysis Buffer (Cell Signalling Technology, Frank- 
furt am M., Germany) supplemented with Protease Inhibitor Mix 
G (Serva, Heidelberg, Germany). After sonication (3 times for 
20 s) and protein quantification (NanoDrop, Thermo Scientific), 
80 |xg/lane of total protein were loaded on a denaturing 10% 
PAA/SDS gel. 

To determine time-dependent basal protein degradation of 
specific NFkB pathway constituents, cells were treated with 
cycloheximide (50|xg/ml). Proteins were isolated at indicated 
time-points. 

For p65: iKB-a co-immunoprecipitation, proteins were isolated 
using NP-40 buffer (50 mM Tris-HCl, 150 mM NaCl, 1% NP- 
40). Three micrograms of p65-specific antibody (clone F6, Santa 
Cruz Biotechnology, Heidelberg, Germany) and 40 |xl G PLUS- 
Agarose (Santa Cruz Biotechnology) were mixed with 1.5 mg 
of total protein and incubated over-night. After centrifugation, 
the pellet was washed three times with HNTG buffer (20 mM 
Hepes pH7.5, 150 mM NaCl, 0.1% Triton X-100, and 10% glyc- 
erol). Subsequently, the pellet was diluted in lx sample buffer 



(100 mM Tris-HCl pH 6.8, 5% glycerol, 0.005% BBP, 5% SDS, 5% 
|3-mercaptoethanol) and incubated at 95°C for 6 min. Pellet was 
pulled down and supernatant was loaded on a 10% PAA/SDS gel. 

After protein transfer, membranes were blocked in lx TBS con- 
taining 0.1% Tween and 5% BSA (Serva, Heidelberg, Germany) or 
5% milk for 30 min. The respective primary antibodies were incu- 
bated over-night at 4°C. After washing with TBST, membranes 
were incubated with specific secondary antibody (1 h at RT) and 
signal detection was performed by ECL (Western Lightning Plus- 
ECL, Perkin Elmer). Signals were digitally documented using the 
Fluorchem FC device (Alpha Innotech/Biozym, Hess. Oldendorf, 
Germany). 

Primary antibodies used in this study were: anti-p65 (clone 
A, dilution for WB: 1:200, Santa Cruz Biotechnology), anti-p65 
(clone F6; dilution for IP: 1:25, Santa Cruz Biotechnology), anti- 
phospho-p65 (Ser536, dilution for WB: 1:500, Santa Cruz Biotech- 
nology), anti-lKB-a (#9242, dilution: 1:750, Cell Signalling Tech- 
nology), anti-phospho-lKBa (Ser32, clone: 14D4, dilution: 1:400, 
Cell Signalling Technology), and anti-actin (clone: C4, dilution: 
1:2000, MP Biomedical, Germany). Secondary antibodies were: 
anti-rabbit IgG HRP-linked antibody (dilution: 1:2000, Cell Sig- 
nalling Technology) and anti-mouse IgG (dilution: 1:2000, Cell 
Signalling Technology). 

2. 1.3. Data acquisition and processing 

Every experiment was replicated at least eight times and raw data 
have been used for modeling purposes. For every protein, four dig- 
ital pictures at different exposure times were taken and signals were 
quantified using the image analysis software Quantity- One (Bio- 
rad, Munich, Germany). Each numerical data point was normal- 
ized to respective actin amounts, which served as loading control 
for the samples. For comparability between different experiments, 
all data points were further normalized to the appropriate control 
on each blot (hepatocytes without TNFa stimulation) to a fixed 
value (c = 100 arbitrary units). 

2.2. COMPUTATIONAL 

A system of ordinary differential equations (ODEs) was set up with 
the software COPASI (Hoops et al., 2006) and integrated using 
the LSODA integrator as implemented in this software (Petzold, 
1983). Sensitivity analysis was employed to determine the influ- 
ence of certain parameters on systems outcomes, e.g., amplitudes 
as described in the text. The sensitivity analysis was also used as 
implemented in COPASI. Parameter fitting was performed using 
the particle swarm algorithm (Kennedy and Eberhart, 1995), again 
as implemented in COPASI. 

3. RESULTS 

3.1. INITIAL TNFcx/NFkB SIGNALING MODEL 

In order to investigate hepatocyte-specific behavior of NFkB sig- 
naling, we developed an initial model based on the most important 
NFkB pathway constituents reported in the literature for hepato- 
cytes (Wullaert et al., 2007; Figure 1A). The involved processes 
were modeled similar to published models (O'Dea et al., 2007; 
Ashall et al, 2009). All processes upstream of the IkB kinase (IKK|3) 
were lumped in the model assuming that they happen very fast 
(Delhase et al., 1999). Phosphorylated IKK|3 is active and phospho- 
rylates free and p65-bound IkBcl Phosphorylated IkBcx dissociates 
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FIGURE 1 | Initial mathematical model of TIMFcx/IMFkB signaling in 
hepatocytes. (A) Schematic representation of the initial model created 
with the software CellDesigner4.2 indicating the considered 
components and reactions (Funahashi et al., 2008). Boxes with round 
corners symbolize proteins, straight boxes indicate DNA, and 
parallelograms represent RNA. (B) Time-course of TN Fa-induced 
signaling and model trajectories. 10 6 primary mouse hepatocytes were 
stimulated with 10ng/mlTNFa for the indicated time. From total cellular 
lysates phosphorylation (pphodpho-serine 536 p65 and phospho-serine 
32 kBa) and total protein (kBa) were examined by immunoblotting. 
Equal loading was confirmed by reprobing of the immunoblots with an 



anti-actin antibody. Quantification was performed using 
chemiluminescence, detection by the Fluorchem FC device and the 
image analysis software Quantity-One. Representative immunoblots 
are indicated above each panel. Quantifications from four independent 
preparations of primary mouse hepatocytes are indicated by symbols. 
The time-course of kBa mRNA induction was determined by real-time 
PCR. Representative results for two independent experiments are 
depicted. See Tables S11 and S12 in Supplementary Material for exact 
values of the experimental data. Black lines indicate model trajectories 
for 30 parameter sets, fitting and simulation was done with the 
software COPASl.The concentrations are given in arbitarty units (a. u.). 



from p65, and is rapidly degraded (Ghosh et al, 1998). Free p65 is 
imported to the nucleus where it induces transcription of different 
target genes such as IkBcx and A20 (Scott et al., 1993; Tian et al., 
2005). In our initial model we introduced stimulated and basal 
transcription, translation, nuclear mRNA export, and degrada- 
tion of IkBoi as well as IicBa-mRNA. Newly synthesized IkBcx can be 
phosphorylated and degraded or it can reenter the nucleus where it 
forms a complex with p65 which is then exported to the cytoplasm 
(Arenzana-Seisdedos et al., 1995). Furthermore, p65 turnover was 
included in form of translation and degradation reactions while 
assuming a constant concentration of the corresponding mRNA. 
Most reactions were modeled as reversible reactions with mass 
action kinetics, while phosphorylation of complexed and free IkBoi 
was modeled by a Michaelis-Menten kinetics with IKK|3 being the 
catalytic enzyme in this term. The included protein synthesis and 
degradation for p65 and IkBcx, and the reversible kinetics as well 
as a basal IKK|3 phosphorylation led to basal levels of all species 
with exception of the ligand TNFa. We assumed the model to be 
in a steady- state before administration of TNFa. For more details 
on the reactions and the parameters see Tables SI, S3, S5-S7 in 
Supplementary Material. 



3.2. PARAMETRIZATION OF THE MODEL BASED ON EXPERIMENTAL 
DATA 

To parametrize the initial model of TNFa-induced NFkB signal- 
ing, quantitative time-courses of NFkB pathway constituents were 
measured after TNFa treatment of primary murine hepatocytes. 
Since it has been shown that TNFa operates very early during 
liver regeneration (Fausto et al., 2006), we investigated the sys- 
tem for 4 h after TNFa stimulation and measured concentrations 
and/or activities using narrow time-point sampling, especially in 
the beginning of this time frame. 

We detected the phosphorylation of iKBa and p65 by 
immunoblotting using total cellular lysates isolated at different 
time-points after TNFa stimulation. Rather constant levels of 
total p65 during the time-course were observed (data not shown). 
However, TNFa stimulation resulted in a rapid increase of phos- 
phorylated p65 after 10 min, followed by a strong decrease until the 
original level was reached after approximately 40 min (Figure IB). 
IKBa was phosphorylated approximately with the same dynamics 
as p65; however, slightly faster and with a second peak of activa- 
tion around 60 min, which lasted until 120 min after stimulation 
(Figure IB). 
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In line with the notion that phosphorylation of IkBcx is 
associated with K48 -linked polyubiquitination and subsequent 
degradation, a rapid reduction of the iKBa protein level from 10 to 
40 min after TNFa treatment and fast recovery after 60 min until 
the end of the experiment were detected (Figure IB). To determine 
to what extent the p65 activity triggered IkBcx gene transcription, 
real-time PCR analysis was performed revealing increased IKBa 
mRNA levels starting from 20 min with a maximum peak level 
around 60 min after stimulation. These experimental data were 
used to perform parameter estimation for the initial model. 

To further constrain the parameters of the model, basal 
turnover of the examined signaling components was measured. 
For this primary mouse hepatocytes were pretreated with the 
inhibitor of protein synthesis cycloheximide (CHX) and time- 
courses up to 360 min without TNFa administration were mea- 
sured for p65 and IkBcx by immunoblotting (Figure 2A). The 
results revealed high p65 stability throughout the entire experi- 
mental time frame, indicating little or no p65 turnover. On the 
other hand, IkBcx showed a 40% degradation within 6h after 
CHX treatment (Figure 2B). Therefore, a constraint was intro- 
duced for parameter estimation to only consider parameter sets 
that ensured a stable level of p65. Individual parameter values 
for the degradation of p65 and IkBcx in their different binding 
states could not be derived from the measurements of total p65 
and IKBa concentrations. Instead, parameters of p65 and IKBa 
degradation were constrained in the parameter fitting procedure 
such that the resulting overall degradation rates were close to the 
measured values. 

The initial model was able to fit most of the data well. The 
fitting, which was computationally very time consuming, was 
repeated 60 times with random start values, the 30 best fits were 



selected for further analysis. However, a more detailed inspec- 
tion of the resulting model trajectories revealed that simultaneous 
fitting of phosphorylated p65, IKBa mRNA, and IkBcx time-courses 
was not sufficiently possible with this preliminary model. More 
specifically, considering the time-courses of phosphorylated p65, 
phosphorylated IKBa, and total IkBcx revealed two distinct clusters 
of parameter sets, none of which could capture all the features of 
the observed data. In the first cluster of parameter sets only the 
initial fast increase of phosphorylated p65 was captured, the fol- 
lowing equally fast decrease however was predicted to be slower 
than the experimental data would suggest; furthermore, IkBcx lev- 
els in the model did not decrease as strong as the experimental data 
would suggest. In the second group of parameter sets the dynamics 
of phosphorylated p65 were captured better, while total IkBcx levels 
went into a too strong overshoot above the steady-state concen- 
tration at approximately 100 min. Also the amplitude of the IKBa 
mRNA dynamics was significantly smaller than the measurements 
suggested. These inconsistencies are due to a tight coupling of free 
p65 and IkBcx levels in the initial model and cannot be resolved 
without changes to the model structure. 

3.3. EXTENDED MODEL INCLUDING ADDITIONAL p65 
PHOSPHORYLATION 

To resolve these discrepancies and improve the description of the 
experimentally observed induction of IkBcx mRNA and subse- 
quent protein levels, we introduced into the initial model reac- 
tions that could modify the transcriptional activity of p65 in the 
nucleus. Since besides phosphorylation at serine 536 additional 
phosphorylation of p65 by nuclear kinases has been proposed, 
which could modify target gene expression including IKBa, we 
considered the activation of a nuclear kinase. The activation of the 
nuclear kinase was modeled generically by three reaction steps with 
mass action kinetics to account for a delay in activation, avoiding 
specific aspects of the upstream activation of the nuclear kinase 
(Figure 3A, Tables S2, S4 in Supplementary Material). 

This extended model was now able to quantitatively fit all exper- 
imental data (Figure 3B). In order to analyze the uniqueness of 
the fit and thus parameter identifiability (Raue et al., 2009), the 
fitting procedure was repeated 30 times with random start values 
resulting in 30 parameter sets fitting the experimental data well. 
Most of the parameters were not identifiable as depicted in Tables 
S8-S10 in Supplementary Material. However, the kinetic constants 
for IKK|3 dephosphorylation and IkBcx basal transcription could 
be narrowed down to a large extent. 

3.4. IMPACT OF PROTEIN TURNOVER ON THE DYNAMICS OF 
TNFk/NFkB SIGNALING 

The role of protein turnover has been neglected in most of the 
studies involving NFkB signaling. However, IkBcx degradation 
might be a critical step because its protein levels immediately 
decreased after activation of the system (IkBcx phosphorylation). 
In addition, the system's main activated effector, p65, is in turn 
able to increase IkBcx protein levels by promoting its transcription, 
thus establishing a negative feedback loop in the system. 

Therefore, we computationally analyzed the impact of protein 
degradation and turnover on important signaling readouts such 
as the amplitude peak of phosphorylated p65 or the integrated 
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FIGURE 2 | Basal protein degradation of NFkB pathway components. 

In the absence of TNFa, primary mouse hepatocytes were treated with 
50|xg/ml cyloheximide and were lyzed at the indicated time-points. (A) 
Representative immunoblot of three independent experiments is shown for 
IkB and p65. (B) Median values and standard deviation of measured p65 
and kBa from three independent experiments. 
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FIGURE 3 | Extended mathematical model of TNFa/NFicB signaling in 
hepatocytes. (A) Schematic representation of the extended model created 
with the software CellDesigner4.2 indicating the additional steps added to the 
initial model. (B) Time-course of TN Fa-induced signaling and model 



trajectories of the extended model. Experimental data for phospho-65, 
phospho-kBa, kBa, and kBa mRNA as shown in Figure 1 is depicted by 
symbols. The predicted model trajectories for 30 parameter sets are indicated 
by black lines. The concentrations are given in arbitarty units (a. u.). 



response (area under curve) which represent the signal strength. 
The ensemble of model parametrizations produced by 30 indepen- 
dent fits was analyzed again. Each of these parametrizations was 
subjected to a sensitivity analysis calculating the impact of spe- 
cific kinetic parameters and protein turnover on the peak-height 
(amplitude) as well as integrated response of the signal (concen- 
tration of NFkB). The results demonstrated that the properties of 
these models were surprisingly robust as depicted in Figure 4 (for 
the respective heat maps of the IkBoi, phosphorylated IkBoi, and 
IkBcx mRNA amplitude and integral see Figures A1-A3). 

These heat maps showed a high consistency within the 30 differ- 
ent fits despite the unidentifiability of most parameters, indicating 
that these sensitivities were valid for most parameter sets and 
therefore most likely for physiological parameter sets. Interest- 
ingly, few parameter sets, e.g., parameter set 30, showed differences 
in comparison with most parameter sets. Analyzing the specific 
behavior of these unique models, revealed a high sensitivity toward 
many parameters, indicating these systems to be in an instable 
state with respect to perturbation. Therefore, it can be assumed 
that these few models, although representing the data well, do not 
correspond to the physiological states in cells. 

When analyzing especially the parameters involved in protein 
synthesis and degradation, which regulate protein level, it became 
apparent that among the prominent parameters influencing the 
integral of the signal were those involved in transferring the TNFa 
stimulation into phosphorylated p65 liberation in form of IKK|3 
phosphorylation (parameters 24, 26, 27 in Figure 4) and p65 
translation and degradation (parameters 42, 37). Interestingly, 
changing translation and degradation of p65 proportionally, so 
that turnover was changed but the protein level stayed constant, 
showed very little impact (parameter 50). Parameters describing 



basal IkBoi transcription, translation, and turnover (parameters 
22, 23, 49) were also influencing the integral drastically, whilst 
its degradation and active transcription (parameters 7-10, 19-21) 
showed very little impact. 

Among the parameters determining the signal amplitude were 
degradation and translation of p65 (parameters 37, 42); again, 
despite its turnover (parameter 49) having very little influence. In 
general, the amplitude of phosphorylated p65 was controlled by 
fewer parameters than the integral of phosphorylated p65. 

With respect to the amplitudes and integrals of phosphorylated 
p65, total iKBa, phosphorylated IkBcx, and iKBa mRNA, the system 
typically depended on the p65 translation and degradation (para- 
meters 37, 42), the activated IKK|3 phosphorylation (parameters 
24, 26), and the IkBcx mRNA degradation in the cytoplasm (para- 
meter 11), whilst changing p65 turnover (assuming unchanged 
steady-state levels of p65, parameter 49) had limited impact. 

3.5. VALIDATION OF THE MODEL AND PREDICTION OF COMPLEX 
FORMATION 

To assess improvements achieved through the development of the 
extended model, we compared the ratios of free IKBa before TNFa 
administration predicted by the initial model and the extended 
model. In the initial model all parameters sets propose that at least 
50% of total IKBa is in an unbound state. On the other hand in the 
extended model the median lies at 20% (Figure 5) which is well 
in line with the experimentally observed 15% reported by Rice 
(1993). 

Finally, we validated our mathematical model of TNFa- 
induced NFkB signaling in primary mouse hepatocytes by using 
experimental data not utilized for parametrization of the model. 
p65:lKBa co-immunoprecipitation was performed to measure 
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i — Parameter sets 

1. It,: complex association 

2. k 2 : complex association 

3. k{. complex nuclear shuttle 

4. k 2 : complex nuclear shuttle 

5. k^: IkBa complex phoshorylation 

6. k M : IkBa complex phoshorylation 

7. k: IkBa degradation complex cytoplasm 

8. k: IkBa degradation complex nucleus 

9. k: IkBa degradation cytoplasm 

10. k: IkBa degradation nucleus 

11. k: IkBa-mRNA degradation cytoplasm 

12. k: IkBa-mRNA degradation nucleus 

13. k: IkBa-mRNA transport 

14. k^ IkBa nuclear shuttle 

15. k 2 : IkBa nuclear shuttle 

16. k: plkBa active degradation 

17. k^: IkBa phosphorylation 

18. k M : IkBa phosphorylation 

19. k A1 : IkBa transcription_pSer 536 p65 

20. k A2 : IkBa transcription_p65 

21. k„ 3 : IkBa transcription_pSer 536 p65 

22. k: IkBa basal transcription 

23. k: IkBa translation 

24. k^: IKKp active phosphorylation 

25. k M : IKKp active phosphorylation 

26. v max : IKKp active phosphorylation 

27. k: IKKp basal phosphorylation 

28. k: IKKp dephosphorylation 

29. k ca ,: Kinase activation 

30. k M : Kinase activation 

31. Kinase deactivation 

32. k: Kinase signal depletion 

33. k: Kinase signal transmission 2 

34. k A : Kinase signal transmission 

35. k M : Kinase signal transmission 

36. v max : Kinase signal transmission 

37. k: p65 degradation 

38. k M : p65 phosphorylation by Kinase 

39. k^: p65 phosphorylation by Kinase 

40. k.,: p65 nuclear shuttle 

41. k 2 : p65 nuclear shuttle 

42. k: p65 translation 

43. k: pSer 536 p65 translation 

44. k^: pp65 phosphorylation by Kinase 

45. k M : pp65 phosphorylation by Kinase 

46. k M : pSer 276 p65 dephosphorylation 

47. TNFa amount 

48. k: TNFa degradation 

49. Turnover IkBa 

50. Turnover p65 



FIGURE 4 | Impact of parameter changes on amplitude and 
integral of phospho-p65. For each of 30 parameter sets the 
sensitivities with respect to every model parameter are shown in a 
heat map. Blue shades denote negative, red shades positive 
sensitivity coefficients. Color shading is scaled from -1.0 to +1.0. The 



parameter sets were clustered hierarchically according to the shown 
sensitivities, however, this is merely to improve visual survey and 
there is no direct correlation to the similarity of parameters. The 
clustering was preformed with the help of the software Cluster 3.0 
(de Hoon et al., 2004). 



complexed p65 concentrations over time (Figure 6B). Since 
most of the parameter values were not identifiable, we cre- 
ated a model ensemble by conducting the parameter fitting 30 
times. All 30 fits, representing different parameter sets all fit- 
ting the experimental data, were used to predict the observed 



experimental co-immunoprecipitation data. For this purpose, the 
kinetic parameters were kept unchanged; only the scaling factors 
were refitted to the validation data (Figure 6A). 

The experimental results showed the initial level of p65:lKBa 
complexes strongly decreased 5 to lOmin after stimulation with 
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FIGURE 5 | Comparison of free IkBk predicted by the initial and the 
extended model. The level of free, unphosphorylated licBafor30 
parameter sets of both models was predicted and is depicted as a box plot. 
According to literature at most 15% of kBa is unbound of p65 (Rice, 1993). 



p65:kBa complex formation 
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FIGURE 6 | Model prediction and experimental validation of p65: licBa 
complex formation. About 10 6 primary mouse hepatocytes were 
stimulated with 10ng/mlTNFa for the indicated time. From NP-40 lysates 
complexes of p65 and kBa were co-immunoprecipitated (IP) using an 
anti-p65 antibody and probing the immunoblots (IB) with an anti-kBa 
antibody. A representative immunoblot is shown and quantifications for two 
independent hepatocyte preparations are indicated by symbols. Black lines 
indicate the time-course predictions from 30 parameter sets. 



TNFa. During this time frame, concentrations of serine 536 
phosphorylated p65 increased and IKBa concentration decreased, 
while after 40 min the I KBa and complex concentrations increased 



again. Between 60 and 120 min the complexed p65:lKBa reached 
their initial level, followed by a second minor decrease. The model 
predicted the full dynamics until the second decrease at approxi- 
mately 180 min, which was likely due to additional feedback loops 
in the system. 

4. DISCUSSION 

The involvement of NFkB signaling is critical in the context 
of liver injury and regeneration as well as inflammation and 
tumorigenesis (Luedde and Schwabe, 2011). Especially during 
chronic liver disease, activation of this pathway triggers hepatocyte 
proliferation and (at later time-points) is eventually supporting 
tumor development and progression. Hepatocellular carcinoma 
(HCC) is the third most common cause of cancer- related death 
worldwide (Breuhahn et al., 2011) and represents a paradigm of 
inflammation-induced cancer, indicating a strong link between 
activation of pro -inflammatory pathways such as TNF-induced 
NFkB signaling and tumorigenesis. However, the impact of NFkB 
signaling on early phases of hepatocarcinogenesis is highly com- 
plex, since this pathway is activated in different cell types including 
hepatocytes (which represent the major cell population), hepatic 
stellate cells, and resident immune cells (Kupffer cells) (Amann 
et al., 2009). In order to understand the impact of this pathway 
in the earliest stages of liver damage and tumorigenesis, it is of 
central relevance to generate mathematical models describing its 
dynamic behavior in normal and non- malignantly transformed 
cells. These tools represent the basis for further comprehensive 
systems biology studies, essential for a deeper understanding of 
liver cancer initiation and progression. 

Our initial model for TNFa-induced NFkB signaling in mouse 
hepatocytes was constructed by applying central aspects of estab- 
lished computational models of NFkB signaling in other organ- 
isms or cell types (O'Dea et al, 2007; Ashall et al, 2009). Our goal 
was a reparametrization of this model using hepatocyte-specific 
experimental data obtained from comprehensive, quantitative, 
and time-resolved protein and transcript analyses. However, the 
model was not able to satisfyingly fit all experimental data simul- 
taneously. Specifically, it was not able to describe the difference 
in the dynamics of phospho-p65 and free IKBa concentrations 
during the first hour after administration of TNFa. While the ini- 
tial rise of phospho-p65 and the corresponding decrease of IKBa 
occur on a similar fast time scale, phospho-p65 returns to near 
its steady-state value faster than IKBa. In the initial model, the 
dephosphorylation of p65 is assumed to coincide with its binding 
to I KBa implying it cannot happen faster than the regeneration 
of IKBa. Since this is contrary to our observations, we uncoupled 
p65 dephosphorylation from its binding to IKBa by introducing 
a more detailed description of p65 phosphorylation and dephos- 
phorylation including a state in which p65 is not phosphorylated 
but also not bound to IKBa in the nucleus. 

It is well known that p65 can be phosphorylated by a number of 
different kinases at different residues such as serine 536, which is 
measured by our immunoblotting approach (Viatour et al., 2005). 
Among the different kinases that may effect serine 536 phosphory- 
lation, only IKK|3 is activated by TNFa-induced signaling (Viatour 
et al., 2005). This is the process that is implicitly included in 
our initial model. The extended model contains in addition the 
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phosphorylation of p65 at serine 276 by a nuclear kinase. Our 
data suggests uncoupling of p65 phosphorylation and its bind- 
ing to IkBoi, which was achieved by implementing an additional 
phosphorylation and corresponding dephosphorylation process 
in the model. Adjusted like this, the model can indeed reproduce 
the observed differences in the dynamics of phospho-p65 and free 
IkBol In addition, the revised model also convincingly reproduced 
the amplitude of the IkBcx mRNA measurements, which was not 
possible using the initial model. 

A potential candidate for such a nuclear kinase could be 
mitogen- and stress-activated protein kinase 1 (MSK1) that has 
also been shown to be expressed in human liver cells (Deak 
et al., 1998). However, MSK1 appears to have no direct impact 
on iKBa expression and the phosphorylation of p65 by MSK1 is 
controversially discussed (Joo and Jetten, 2008). 

The extended model is able to explain features of NFkB signal- 
ing in hepatocytes as displayed by our measurements. However, 
NFkB signaling is known to involve more constituents and regula- 
tions than considered in this study. Therefore, other components 
can be included in future research. Tumor necrosis factor a- 
induced protein 3 (TNFAIP3/A20) is one considerable protein, 
involved in a negative loop inhibiting p65 phosphorylation acting 
upstream of IKK-|3 (Opipari et al., 1990). This protein, together 
with IkBcx is probably the most important factor in the acquisition 
of phosphorylated p65 oscillatory behavior critical for transcrip- 
tional target fate (Ashall et al., 2009). Another NFkB constituent 
involved in negative feedback loop with high impact in NFkB sig- 
naling oscillations is IkBs. This protein belongs to the IkB inhibitor 
family of signaling whose activity has been detected later and out 
of phase if compared with IkBcx. This special feature leads IkBs 
to dampen the effect of IkBcx during the last phase of sustained 
NFkB activity and to mediate the terminations of the signaling in 
response to transient stimulation (Kearns et al., 2006). From our 
preliminary data, both protein TNFAIP3/A20 and IkBs present 
a significant expression in primary hepatocytes following TNFa 
stimulation and this data will be used for parametrization of fur- 
ther improved models. Moreover signals coming not directly from 
TNFa but able to fine-tune the system especially in the nucleus 
should be taken into account. 

A vital methodological question remains in the context of this 
work: What is the justification to claim that the model proposed 
here is a model for NFkB signaling in mouse hepatocytes? We have 
consulted literature to make sure all the processes and compo- 
nents included in the model are actually present in mammalian 
hepatocytes. On the other hand, to our knowledge none of these 
components or processes is actually exclusive to liver cells or hepa- 
tocytes. Basically this leaves the parametrization of the model as 
being hepatocyte-specific. Parameter estimation was performed 
with data measured in murine hepatocytes and the presented 
model was able to fit the data satisfyingly. However, some of 
the parameters are unidentifiable with the available data, i.e., for 
different values of the parameter the model fits the data equally 
well. Instead of eliminating the unidentifiabilities - by performing 
additional experiments and/or by simplifying the model - here 
we focused on making predictions from the model that are valid 
for all parametrizations of the model that are compatible with 
the experimental data. As a pragmatic way to characterize the 



set of parameter values that let the model fit the data, a heuris- 
tic sampling of parameter space was performed, a technique that 
was recently applied in other studies (Maiwald et al., 2010; Lev- 
ering et al., 2012; Wegner et al., 2012). The parameter estimation 
process was repeated 30 times, each time starting with different 
random parameter values as initial guess, resulting in 30 differ- 
ent sets of parameters that are compatible with the experimental 
data. This immediately results in a good estimate for the iden- 
tifiability of the different parameters: the more the parameters 
differ between the repeated parameter estimation runs the less 
identifiable they are. 

For drawing conclusions from our model or creating predic- 
tions for validation we need to make sure that these conclusions 
or predictions are valid for all or at least most of the 30 parameter 
sets. Indeed, the model validation described in the Results section 
proved to be conserved for all 30 parameter sets. Similarly the sen- 
sitivities that were calculated in order to study the role of protein 
turnover turned out to be consistent for 25 of the 30 parameter 
sets. While of course we cannot be absolutely certain from using 
only a limited number of samples, we can confidently claim that 
typically the model will predict a certain behavior of NFkB signal- 
ing in hepatocytes whenever the parameters are chosen so that it 
fits the data. 

While our model shares its basic structure with most relevant 
NFkB models in literature, it is based on NFkB dynamics fol- 
lowing TNFa stimulation in murine primary hepatocytes. Our 
experimental setup allows us to disturb hepatocytes with differ- 
ent combination of ligands and to compare data with prediction 
from our model. Different models for different signaling path- 
ways essential for liver homeostasis and response to injury are 
currently being developed, e.g., within the German Virtual Liver 
Network (VLN). In this context our model can be a useful tool for 
comparison and integration of different models. Thus, we have 
established a mathematical model for TNFa/NFKB signaling in 
primary mouse hepatocytes that indicated the importance of basal 
turnover of the pathway components for the dynamic behavior 
and facilitated the prediction of the kinetics of complex formation 
in the pathway. Taken together the model provides an important 
building block to unravel regulation of liver regeneration and the 
impact of inflammatory responses across the scales. 
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Table S1 | Ordinary differential equations describing all species in the 
original model. As initial concentrations the steady-state concentrations were 
used. 

Table S2 | Ordinary differential equations describing all species in the 
model expanded by p65 phosphorylaiton. 
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Table S3 | Parameters and fitting boundaries of the parameters of the initial 
model. 

Table S4 | Parameters and fitting boundaries of the parameters of the 
expanded model. 

Table S5 | Objective function value and exact values of all parameters fitted 
via particle swarm in the parameter estimation starting from random start 
values for the initial model part I. 

Table S6 | Objective function value and exact values of all parameters fitted 
via particle swarm in the parameter estimation starting from random start 
values for the initial model part II. 

Table S7 | Objective function value and exact values of all parameters fitted 
via particle swarm in the parameter estimation starting from random start 
values for the initial model part III. 



Table SB | Objective function value and exact values of all parameters fitted 
via particle swarm in the parameter estimation starting from random start 
values for the expanded model part I. 

Table S9 | Objective function value and exact values of all parameters fitted 
via particle swarm in the parameter estimation starting from random start 
values for the expanded model part II. 

Table S10 | Objective function value and exact values of all parameters 
fitted via particle swarm in the parameter estimation starting from 
random start values for the expanded model part III. 

Table S11 | Relative protein concentrations over time determined from 
immuno blotting. 

Table S12 | Relative mRNA concentrations over time determined from 
quantitative reverse transcription polymerase chain reaction. 



REFERENCES 

Amann, T., Bataille, E, Spruss, T., 
Miihlbauer, M., Gabele, E., 
Scholmerich, J., et al. (2009). 
Activated hepatic stellate cells 
promote tumorigenicity of hepato- 
cellular carcinoma. Cancer Sci. 100, 
646-653. 

Arenzana-Seisdedos, E, Thompson, J., 
Rodriguez, M. S., Bachelerie, E, 
Thoma, D., and Hay, R. T. (1995). 
Inducible nuclear expression of 
newly synthesized Ik Ba negatively 
regulates dna-binding and tran- 
scriptional activities of NFk-B. Mol. 
Cell. Biol. 15, 2689-2696. 

Ashall, L., Horton, C. A., Nelson, D. E., 
Paszek, P., Harper, C. V., Sillitoe, K., 
et al. (2009). Pulsatile stimulation 
determines timing and specificity 
of NFK-B-dependent transcription. 
Science 324, 242-246. 

Bohuslav, J., Chen, L.-E, Kwon, H., 
Mu, Y., and Greene, W. C. (2004). 
p53 induces NFk-B activation by an 
Ik B kinase-independent mechanism 
involving phosphorylation of p65 by 
ribosomal S6 kinase 1 . /. Biol. Chem. 
279, 26115-26125. 

Breuhahn, K., Gores, G., and Schir- 
macher, P. (2011). Strategies for 
hepatocellular carcinoma therapy 
and diagnostics: lessons learned 
from high throughput and pro- 
filing approaches. Hepatology 53, 
2112-2121. 

Carlotti, E, Chapman, R., Dower, S. 
K., and Qwarnstrom, E. E. (1999). 
Activation of nuclear factor k (B 
in single living cells, dependence 
of nuclear translocation and anti- 
apoptotic function on egfprela con- 
centration. /. Biol. Chem. 274, 
37941-37949. 

Carlotti, E, Dower, S. K., and Qwarn- 
strom, E. E. (2000). Dynamic shut- 
tling of nuclear factor kB between 
the nucleus and cytoplasm as a con- 
sequence of inhibitor dissociation. /. 
Biol. Chem. 275,41028-41034. 



Castoldi, M., Vujic Spasic, M., Alta- 
mura, S., Elmen, J., Lindow, M., Kiss, 
J., et al. (2011). The liver-specific 
microrna mir-122 controls systemic 
iron homeostasis in mice. /. Clin. 
Invest. 121, 1386-1396. 

Cheong, R., Bergmann, A., Werner, S. 
L., Regal, J., Hoffmann, A., and 
Levchenko, A. (2006). Transient IkB 
kinase activity mediates temporal 
NFk-B dynamics in response to 
a wide range of tumor necrosis 
factora- doses. /. Biol. Chem. 281, 
2945-2950. 

de Hoon, M., Imoto, S., Nolan, J., and 
Miyano, S. (2004). Open source clus- 
tering software. Bioinformatics 20, 
1453-1454. 

Deak, M., Clifton, A. D., Lucocq, L. M., 
and Alessi, D. R. (1998). Mitogen- 
and stress -activated protein kinase- 
1 (MSK1) is directly activated by 
MAPK and SAPK2/p38, and may 
mediate activation of CREB. EMBO 
}. 17,4426-4441. 

Delhase, M., Hayakawa, M., Chen, Y., 
and Karin, M. (1999). Positive and 
negative regulation of IkB kinase 
activity through IKKp subunit phos- 
phorylation. Science 284, 309-313. 

Diehl, A. M., and Rai, R. (1996). Review: 
regulation of liver regeneration by 
pro -inflammatory cytokines. /. Gas- 
troenterol. Hepatol. 11,466-470. 

Duran, A., Diaz-Meco, M. T, and 
Moscat, J. (2003). Essential role 
of RelA Ser3 1 1 phosphorylation by 
yPKC in NFk-B transcriptional acti- 
vation. EMBO }. 22, 3910-3918. 

Fausto, N., Campbell, J. S., and Riehle, K. 
J. (2006). Liver regeneration. Hepa- 
tology 43(Suppl. 1), S45-S53. 

Funahashi, A., Matsuoka, Y., Jouraku, 

A. , Morohashi, M., Kikuchi, N., and 
Kitano, H. (2008). CellDesigner 3.5: 
a versatile modeling tool for bio- 
chemical networks. Proc. IEEE 96, 
1254-1265. 

Ghosh, S., May, M. J., and Kopp, E. 

B. (1998). NFk-B and Rel proteins: 



evolutionarily conserved mediators 
of immune responses. Annu. Rev. 
Immunol. 16, 225-260. 
Hayden, M. S., and Ghosh, S. (2012). 
NFk-B, the first quarter-century: 
remarkable progress and outstand- 
ing questions. Genes Dev. 26, 
203-234. 

Hoffmann, A., Levchenko, A., Scott, M. 
L., and Baltimore, D. (2002). The Ik 
B-NFk-B signaling module: tempo- 
ral control and selective gene activa- 
tion. Science 298, 1241-1245. 

Hoops, S., Sahle, S., Gauges, R., Lee, C, 
Pahle, J., Simus, N., et al. (2006). 
COPASI-a complex pathway simu- 
lator. Bioinformatics 22, 3067-3074. 

Ihekwaba, A. E. C, Broomhead, D. 
S., Grimley, R. L., Benson, N., and 
Kell, D. B. (2004). Sensitivity analysis 
of parameters controlling oscillatory 
signalling in the NFk-B pathway: the 
roles of ikk and Ik Ba. Syst. Biol. 
(Stevenage) 1, 93-103. 

Jamaluddin, M., Wang, S., Boldogh, 
I., Tian, B., and Brasier, A. R. 
(2007). TNFot-induced NFk- 
B/RelA Ser(276) phosphorylation 
and enhanceosome formation is 
mediated by an ROS-dependent 
PKAc pathway. Cell. Signal. 19, 
1419-1433. 

Joo, J. H., and Jetten, A. M. (2008). NFK- 
B-dependent transcriptional activa- 
tion in lung carcinoma cells by 
farnesol involves p65/RelA(Ser276) 
phosphorylation via the MEK- 
MSK1 signaling pathway. /. Biol. 
Chem. 283,16391-16399. 

Kearns, J. D., Basak, S., Werner, S. L., 
Huang, C. S., and Hoffmann, A. 
(2006). IkBoi provides negative feed- 
back to control NFk-B oscillations, 
signaling dynamics, and inflamma- 
tory gene expression. /. CellBiol. 173, 
659-664. 

Kennedy, J., and Eberhart, R. (1995). 
Particle swarm optimization. Proc. 
Int. It. Conf. Neural Netw. 4, 
1942-1948. 



Krappmann, D., Wulczyn, F. G., and 
Scheidereit, C. (1996). Different 
mechanisms control signal-induced 
degradation and basal turnover of 
the NF-kappaB inhibitor IkB alpha 
in vivo. EMBO /. 15, 6716-6726. 

Levering, J., Musters, M. W. J. M., 
Bekker, M., Bellomo, D., Fiedler, T, 
de Vos, W. M., et al. (2012). Role of 
phosphate in the central metabolism 
of two lactic acid bacteria — a com- 
parative systems biology approach. 
FEBS I. 279, 1274-1290. 

Li, J., Campbell, J. S., Mitchell, C, 
McMahan, R. S., Yu, X., Riehle, K. J., 
et al. (2009). Relationships between 
deficits in tissue mass and transcrip- 
tional programs after partial hepa- 
tectomy in mice. Am. /. Pathol. 175, 
947-957. 

Luedde, T, and Schwabe, R. F. (2011). 
NFk-B in the liver - linking injury, 
fibrosis and hepatocellular carci- 
noma. Nat. Rev. Gastroenterol. Hepa- 
tol. 8, 108-118. 

Madrid, L. V., Mayo, M. W, Reuther, 
J. Y, and Baldwin, A. S. (2001). 
Akt stimulates the transactivation 
potential of the RelA/p65 Subunit 
of NF-kappa B through utilization 
of the Ik (B kinase and activa- 
tion of the mitogen-activated pro- 
tein kinase p38. /. Biol. Chem. 276, 
18934-18940. 

Maiwald, T, Schneider, A., Busch, H., 
Sahle, S., Gretz, N., Weiss, T. S., 
et al. (2010). Combining theoreti- 
cal analysis and experimental data 
generation reveals IRF9 as a cru- 
cial factor for accelerating interferon 
a -induced early antiviral signalling. 
FEBS J. 277, 4741-4754. 

Mathes, E., O'Dea, E. L., Hoffmann, 
A., and Ghosh, G. (2008). NFk-B 
dictates the degradation pathway of 
IkBoi. EMBO J. 27, 1357-1367. 

Michalopoulos, G. K. (2010). Liver 
regeneration after partial hepatec- 
tomy: critical analysis of mechanistic 
dilemmas. Am. J. Pathol. 176, 2-13. 



Frontiers in Physiology | Systems Biology 



December 2012 | Volume 3 | Article 466 | 10 



Pinna et al. 



NFkB signaling in primary hepatocytes 



Nelson, D. E., Ihekwaba, A. E. C, Elliott, 
M., Johnson, J. R., Gibney, C. A., 
Foreman, B. E., et al. (2004). Oscilla- 
tions in NFk-B signaling control the 
dynamics of gene expression. Science 
306, 704-708. 

O'Dea, E. L., Barken, D., Peralta, R. Q., 
Tran, K. T., Werner, S. L., Kearns, 
J. D., et al. (2007). A homeostatic 
model of Ik B metabolism to con- 
trol constitutive NFk-B activity. Mol. 
Syst.Biol. 3,111. 

Opipari, A. W., Boguski, M. S., and 
Dixit, V. M. (1990). The A20 cDNA 
induced by tumor necrosis factor 
a encodes a novel type of zinc 
finger protein. /. Biol. Chem. 265, 
14705-14708. 

Pando, M. P., and Verma, I. M. (2000). 
Signal-dependent and -independent 
degradation of free and NFk-B- 
bound Ik Ba. /. Biol. Chem. 275, 
21278-21286. 

Papp, V., Dezso, K., Laszlo, V., Nagy, 
P., and Paku, S. (2009). Archi- 
tectural changes during regenera- 
tive and ontogenic liver growth 
in the rat. Liver Transpl. 15, 
177-183. 

Petzold, L. (1983). Automatic selection 
of methods for solving stiff and non- 
stiff systems of ordinary differential 
equations. SIAMJ. Sci. Stat. Comput. 
4, 136-148. 

Raue, A., Kreutz, C, Maiwald, T., Bach- 
mann, J., Schilling, M., Klingmiiller, 
U., et al. (2009). Structural and prac- 
tical identifiability analysis of par- 
tially observed dynamical models 



by exploiting the profile likelihood. 
Bioinformatics 25, 1923-1929. 

Rice, N. R. and Ernst, M. K. (1993). 
In vivo control of NFk-B activation 
by IkBoi. EMBO }. 12, 4685-4695. 

Sakurai, H., Chiba, H., Miyoshi, 
H., Sugita, T., and Toriumi, W. 
(1999). IkB kinases phospho- 
rylate NFk-B p65 subunit on 
serine 536 in the transactiva- 
tion domain. /. Biol. Chem. 274, 
30353-30356. 

Scott, M. L., Fujita, T., Liou, H.-C, 
Nolan, G. P., and Baltimore, D. 
(1993). The p65 subunit of NFk-B 
regulates IkB by two distinct mech- 
anisms. Genes Dev. 7, 1266-1276. 

Spooren, A., Kolmus, K., Vermeulen, 
L., Van Wesemael, K., Haegeman, 
G., and Gerlo, S. (2010). Hunt- 
ing for serine 276-phosphorylated 
p65. /. Biomed. Biotechnol. 2010, 
275892. 

Sung, M.-H., Salvatore, L., De Lorenzi, 
R., Indrawan, A., Pasparakis, M., 
Hager, G. L., et al. (2009). Sustained 
oscillations of NFk-B produce dis- 
tinct genome scanning and gene 
expression profiles. PLoS ONE 
4:e7163. doi:10.1371/journal.pone. 
0007163 

Tian, B., Nowak, D. E., Jamalud- 
din, M., Wang, S., and Brasier, A. 
R. (2005). Identification of direct 
genomic targets downstream of 
the nuclear factorK-B transcription 
factor mediating tumor necrosis 
factor signaling. /. Biol. Chem. 280, 
17435-17448. 



Vermeulen, L., De Wilde, G., Van 
Damme, P., Vanden Berghe, W, and 
Haegeman, G. (2003). Transcrip- 
tional activation of the NFk-B p65 
subunit by mitogen- and stress- 
activated protein kinase- 1 (MSK1). 
EMBO J. 22, 1313-1324. 

Viatour, P., Merville, M.-R, Bours, V, 
and Chariot, A. (2005). Phosphory- 
lation of NFk-B and IkB proteins: 
implications in cancer and inflam- 
mation. Trends Biochem. Sci. 30, 
43-52. 

Wang, Y., Paszek, P., Horton, C., 
and Kell, D. (2011). Interactions 
among oscillatory pathways in NFk- 
B signalling. BMC Syst. Biol. 5:23. 
doi:10.1186/1752-0509-5-23 

Wegner, K., Bachmann, A., Schad, J.-U., 
Lucarelli, P., Sahle, S., Nickel, P., et 
al. (2012). Dynamics and feedback 
loops in the transforming growth 
factorp signaling pathway. Biophys. 
Chem. 162, 22-34. 

Werner, S. L., Kearns, J. D., Zadorozh- 
naya, V, Lynch, C., O'Dea, E., Boldin, 
M. P., et al. (2008). Encoding NFk-B 
temporal control in response to tnf: 
distinct roles for the negative regu- 
lators IkBoi and a20. Genes Dev. 22, 
2093-2101. 

Wullaert, A., van Loo, G., Heyninck, 
K., and Beyaert, R. (2007). Hepatic 
tumor necrosis factor signaling and 
nuclear factorK-b: effects on liver 
homeostasis and beyond. Endocr. 
Rev. 28, 365-386. 

Yamada, Y., Kirillova, I., Peschon, J. J., 
and Fausto, N. (1997). Initiation of 



liver growth by tumor necrosis fac- 
tor: deficient liver regeneration in 
mice lacking type i tumor necrosis 
factor receptor. Proc. Natl. Acad. Sci. 
U.S.A. 94, 1441-1446. 

Conflict of Interest Statement: The 

authors declare that the research was 
conducted in the absence of any com- 
mercial or financial relationships that 
could be construed as a potential con- 
flict of interest. 

Received: 11 August 2012; accepted: 26 
November 2012; published online: 31 
December 2012. 

Citation: Pinna F, Sahle S, Beuke K, 
Bissinger M, Tuncay S, D'Alessandro LA, 
Gauges R, Raue A, Timmer J, Kling- 
miiller U, Schirmacher P, Kummer U 
and Breuhahn K (2012) A systems biol- 
ogy study on NFkB signaling in primary 
mouse hepatocytes. Front. Physio. 3:466. 
doi: 10.3389/fphys.2012.00466 
This article was submitted to Frontiers in 
Systems Biology, a specialty of Frontiers 
in Physiology. 

Copyright © 2012 Pinna, Sahle, Beuke, 
Bissinger, Tuncay, D'Alessandro, Gauges, 
Raue, Timmer, Klingmiiller, Schirma- 
cher, Kummer and Breuhahn. This is an 
open-access article distributed under the 
terms of the Creative Commons Attribu- 
tion License, which permits use, distrib- 
ution and reproduction in other forums, 
provided the original authors and source 
are credited and subject to any copy- 
right notices concerning any third-party 
graphics etc. 



Research Topic: 

From structural to molecular systems biology: experimental and 
computational approaches to unravel mechanisms of kinase activity 
regulation in cancer and neurodegeneration 



Topic Editor: 
Matteo Barberis, 

Humboldt University Berlin, Germany; 

Max Planck Institute for Molecular Genetics, Berlin, Germany 



www.frontiersin.org 



December 2012 | Volume 3 | Article 466 | 11 



Pinna et al. 



NFkB signaling in primary hepatocytes 



APPENDIX 



IkBq Amplitude 



IkBq Integral 



1 2 3 4 5 6 7 8 9 1 0 11 12 13 1415 16 17 18 19 2021 22 2324 2526 2728 2930 1 2 3 4 5 6 7 8 9 10 11 12 13 1415 16 17 18 19 2021 22 2324 2526 2728 2930 



Iri 



■ 



I 




< — Parameter sets 

1. complex association 

2. k 2 : complex association 

3. complex nuclear shuttle 

4. k 2 : complex nuclear shuttle 

5. k caf : IkBa complex phoshorylation 

6. k M : IkBa complex phoshorylation 

7. k: IkBa degradation complex cytoplasm 

8. k: IkBa degradation complex nucleus 

9. k: IkBa degradation cytoplasm 

10. k: IkBa degradation nucleus 

11. k: IkBa -mRNA degradation cytoplasm 

12. k: IkBa-mRNA degradation nucleus 

13. k: IkBa-mRNA transport 

14. k.,: IkBa nuclear shuttle 

15. k 2 : IkBa nuclear shuttle 

16. k: plkBa active degradation 

17. k caf : IkBa phosphorylation 

18. k M : IkBa phosphorylation 

19. k^: IkBa transcription_pSer 536 p65 

20. k^: IkBa transcription_p65 

21. k A3 : IkBa transcription_pSer 536 p65 

22. k: IkBa basal transcription 

23. k: IkBa translation 

24. k^: IKKp active phosphorylation 

25. k M : IKKp active phosphorylation 

26. v max : IKKp active phosphorylation 

27. k: IKKp basal phosphorylation 

28. k: IKKp dephosphorylation 

29. k caf : Kinase activation 

30. k M : Kinase activation 

31. Kinase deactivation 

32. k: Kinase signal depletion 

33. k: Kinase signal transmission 2 

34. k A : Kinase signal transmission 
35. 
36. 

37. k: p65 degradation 

38. k M : p65 phosphorylation by Kinase 

39. k caf : p65 phosphorylation by Kinase 

40. k^ p65 nuclear shuttle 
41. 
42. 

43. k: pSer 536 p65 translation 

44. k caf : pp65 phosphorylation by Kinase 

45. k M : pp65 phosphorylation by Kinase 

46. k M : pSer 276 p65 dephosphorylation 

47. TNFa amount 

48. k: TNFa degradation 

49. Turnover IkBa 

50. Turnover p65 



k M : Kinase signal transmission 
v max : Kinase signal transmission 



k 2 : p65 nuclear shuttle 
k: p65 translation 



FIGURE A1 | Impact of parameter changes on amplitude and integral of 
IkBcx. For each of 30 parameter sets the sensitivities with respect to every 
model parameter are shown in a heat map. Blue shades denote negative, red 
shades positive sensitivity coefficients. Color shading is scaled from -1.0 to 



+ 1.0. The parameter sets were clustered hierarchically according to the 
shown sensitivities, however, this is merely to improve visual survay and 
there is no direct correlation to the similarity of parameters. The clustering 
was preformed with the help of the software Cluster 3.0 (Ashall et al., 2009). 
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FIGURE A2 | Impact of parameter changes on amplitude and integral of 
phospho-licBa. For each of 30 parameter sets the sensitivities with respect 
to every model parameter are shown in a heat map. Blue shades denote 
negative, red shades positive sensitivity coefficients. Color shading is scaled 



from -1.0 to +1.0. The parameter sets were clustered hierarchically according 
to the shown sensitivities, however, this is merely to improve visual survey 
and there is no direct correlation to the similarity of parameters. The clustering 
was preformed with the help of the software Cluster 3.0 (Ashall et al., 2009). 
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FIGURE A3 | Impact of parameter changes on amplitude and integral of 
IkBo; mRNA. For each of 30 parameter sets the sensitivities with respect to 
every model parameter are shown in a heat map. Blue shades denote 
negative, red shades positive sensitivity coefficients. Color shading is scaled 



from -1.0 to +1.0. The parameter sets were clustered hierarchically according 
to the shown sensitivities, however, this is merely to improve visual survey 
and there is no direct correlation to the similarity of parameters. The clustering 
was preformed with the help of the software Cluster 3.0 (Ashall et al., 2009). 
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